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1. Introduction 

In supersymmetryljl]], /^-parity violation is one of the possible scenario in beyond standard 
model physics. With many interest, we consider only on the k!j k couplings which arises in lepton 
number violating term in /?-parity violation superpotential given below 

W = pLiLiH 2 + XijkLiLjE c h + X! jk LiQjD c k + X^UfD'pl , (1.1) 

where L,- and Qj are the SU (2)-doublet lepton and quark superfields, Ef ,U[ ,D- the singlet super- 
fields and Hi the Higgs superfields. The subscripts i,j,k are generational indices. Note that 
is antisymmetric under the interchange of the first two indices and X"- k is antisymmetric under the 
interchange of the last two. The first three terms in eqn.(|l.l[) violate lepton number (L) and the last 
term violates baryon number (B) conservation. 

Recently ATLAS group have been studied resonant production of heavy neutral scalar like 
sneutrino and subsequent decay to e\l final state. In their analysis, they put the bounds on sneutrino 
masses (see Ref.|Q]) on the basis of leading order (LO) result. In tevatron, both CDF[^] and D0[^] 
collaboration analyse their data (Run-I as well as Run-II data) using our first results [|5|] on the next 
to leading order (NLO) QCD corrections to sneutrino and charged slepton productions at hadron 
colliders. In their analysis to set bound on these R-parity violating couplings, cross section for 
SM background processes namely Drell-Yan production of pair of leptons (say l + l~ , rv) (see first 



two papers of fll6Q) was considered at the next to next to leading order (NNLO) level while for the 
R-parity violating effects only NLO corrected cross section was used. It was found that the NLO 
QCD effects were quite large ~ 10% — 40% at both Tevatron as well as LHC Therefore, it is de- 
sirable to compute the cross sections for the resonant sneutrino and/or charged slepton productions 
at NNLO in QCD. These results will quantitatively improve the analysis based on high statistics 
data available in the ongoing and future experiments. From the theoretical point of view, higher 
order radiative corrections provide a test of the convergence of the perturbation theory and hence 
the reliable comparison of data with the theory predictions is possible. The fixed order perturbative 
results most often suffer from large uncertainties due to the presence of renormalisation and factori- 
sation scales. They get reduced as we include more and more terms in the perturbative expansion 
thanks to renormalisation group invariance. In this article we have systematically included its scale 
dependence through the renormalisation group equations and we discussed the impact of it in the 
next sections. 

2. Brief discussion of NNLO calculations 

In this section, we describe very briefly, the computation of second order (a s 2 ) QCD radiative 
corrections to resonant production, in hadron colliders, of a sneutrino/charged slepton. We present 
our results in such a way that they can be used for any scalar-pseudoscalar production which is the 
main goal of this work. The inclusive hadronic cross section for the reaction 

H 1 (P l )+H 2 (P 2 )^^(p^)+X, (2.1) 

is given by 

nX' 2 (lll) ^ f 1 dx\ f 1 dx 2 , . 2 w, 2na ( T 2 2 2 
l-^O „u s Jt X\ Jt/Xj X 2 V-^l^ 



a,b=q,q,g 
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with T = ^ , S = (P 1 +P 2 ) 2 , p\=m\, (2.2) 

where H\ and H2 denote the incoming hadrons and X represents an inclusive hadronic state. The 
parton densities denoted by f c (xi,pLp) (c = q,q,g) depend on the scaling variables x\ (i =1,2) 
through pi =XiPi and the mass factorization scale jXp. Here p { (i = 1,2) are the momenta of incom- 
ing partons namely quarks, antiquarks and gluons. The coupling constant X' gets renormalised at 
the renormalisation scale }Xr due to ultraviolet singularities present in the theory. The factorisation 
scale is introduced on the right hand side of the above equation to separate long distant dynamics 
from the perturbatively calculable short distant partonic coefficient functions A ab . A ab depends on 
both pLR and jip in such a way that the entire scale dependence goes away to all orders in pertur- 
bation theory when convoluted with appropriate parton densities. This is due to the fact that the 
observable on the left hand side of the above equation is renormalisation group (RG) invariant with 
respect to both the scales. This implies 



M 2 -j^r = 0, Ai = MF,MR, (2-3) 



Ml 



dps 
d 



X' 2 (pj)A ab {x,m^,f4,f4) 



0. (2.4) 



The partonic coefficient functions that appear in eqn.( |2.2| ) are computable in perturbative QCD 
in terms of strong coupling constant g s . The ultraviolet singularities present in the theory are 
regularised in dimensional regularisation and are removed in MS scheme, introducing the renor- 
malisation scale }JLr at every order in perturbative expansion. In addition, the Yukawa coupling X' 
also gets renormalised due to strong interaction dynamics. Hence, for our computation, we require 
only two renormalisation constants to obtain UV finite partonic coefficient functions, A ab . These 
constants are denoted by Z(/i#) and Z%i(Hr), where the former renormalises the strong coupling 
constant g s and the later Yukawa coupling X' and both the couplings a s (= g s / (4k)) (and X') evolve 
with scale to NNLO through renormalisation group equations: 

d 00 

Hi— jhyX'(nl) = -jy,0*l) Yi-i- (2-5) 

where coefficients j3,- for i = 0, 3 can be found in ^ for SU (N) QCD. The anomalous dimensions 
Yi for i = 0, ...,3 can be obtained from the quark mass anomalous dimensions given in |{7|]. The 
perturbatively calculable A a [, can be expanded in powers of strong coupling constant a s (Hg) as 



i=0 



Aab g ets contributions from various partonic reactions. 

The calculation of various contributions from the partonic reactions involves careful handling 
of divergences that result from one[|8|] and two loop[]|] integrations in the virtual processes and 
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two and three body phase space integrations in the real emission processes. The loop integrals 
often give ultraviolet, soft and collinear divergences. But the phase space integrals give only soft 
and collinear singularities. Soft divergences arise when the momenta of the gluons become zero 
while the collinear diverges arise due to the presence of massless partons. We have regulated all 
the integrals in dimensional regularisation with space time dimension n = 4 + e. The singularities 
manifest themselves as poles in £. 

We have reduced all the one loop tensorial integrals to scalar integrals using the method of 
Passarino-Veltman [ |To| ] in 4 + s dimensions and evaluated resultant scalar integrals exactly. The 
2-loop form factor, ^(ml,/! 2 ), is calculated using the dispersion technique [11]. Two and three 
body phase space integrals are done by choosing appropriate Lorentz frames[12]. Since we inte- 
grate over the total phase space the integrals are Lorentz invariant and therefore frame independent. 
Several routines are made using the algebraic manipulation program FORM[|l3|] in order to perform 
tensorial reduction of one loop integrals and two and three body phase space integrals. 

The UV singularities go away after performing renormalisation through the constants Z and 
Z^i. The soft singularities cancel among virtual and real emission processes [[14]] at every order in 
perturbation theory. The remaining collinear singularities are renormalised systematically using 
mass factorisation[jl5|]. For more details on the computation of NNLO QCD corrections to pro- 



cess of the kind considered here can be found in [16, 17]. The full analytical results for NNLO 
calculation for sneutrino and/or charge slepton can be found out in our original paperrtTj 



3. Results and Discussion 

We considered only the contributions from the first generation of quarks. Since at hadron 
colliders, the resonant production is through the interaction term X[- k LiQjD\ in the Lagrangian 
(see eq.(pTT|)), for j,k = 2,3, the production rate will be suppressed due to the low flux of the 
sea quarks. To obtain the production cross section to a particular order, one has to convolute the 
partonic coefficient functions A a b with the corresponding partem densities f a , both to the same order. 
Further the coupling constants ^ (/Ir) and should also be evaluated using the corresponding 

RGEs (eqn.(2.5)) computed to the same order (more details see Ref[[7|, 18] ). We have used the latest 



MSTW partem densities [[19(] in our numerical code and the corresponding values of a s {Mz) for LO, 
NLO and NNLO provided with the sets. Since we are considering one A/ n non-zero, the LO and 
NLO cross sections get contributions only from dd , dg and dg initiated subprocesses and no other 
quark (antiquark) flavors contribute to this order. At NNLO level, the incoming quarks other than d 
type quarks can also contribute. The total sneutrino production cross section as function of its mass 
is plotted in fig. |] for LHC (left panel) and Run II of Tevatron (right panel) energies. We have set 
the renormalisation scale to be the mass of the sneutrino, pLR = m v . The pair of lines corresponds to 
the two extreme choices of factorisation scale: pLp = 10 m v (upper) and }ip = m v / 10 (lower). The 
plots clearly demonstrate that the NNLO contributions reduce the factorisation scale dependence 
improving the theoretical predictions for sneutrino production cross section. 

The cross section falls off with the sneutrino mass due to the availability of phase space with 
respect to the mass, the choice of /i# = m v and the parton densities. The latter effect, understand- 
ably, is more pronounced at the Tevatron than at the LHC. 
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Figure 1: Total cross-section for the V production as a function ofm v . For smaller values of sneutrino mass 
the upper (lower) set of lines correspond to the factorisation scale )lf = lOmy(O.lmy). For larger values of 
sneutrino mass the lines cross each other. 



In order to estimate the magnitude of the QCD corrections at NLO and NNLO, we define the 
K-factors as follows: 



K W _ a t I -A K (2) _ / 

A — u fof,NLO/ u tot, LO ' A — u tot, NNLO / u tot,LO ■ 

In fig.||, we have plotted both (i = 1,2) as a function of sneutrino mass. We have chosen 
jXf = Hr = niy for this study. At the LHC, The K^ 1 ' varies between 1.23 to 1.46 and between 
1.27 to 1.52 in the mass range 100 GeV < m v < 1 TeV. At the Tevatron, we find that varies 
between 1.55 to 1.53 and between 1.65 to 1.85 for the same mass range. Note that numbers 
for differ from those given in our earlier work [Q] due to the running of A' in the present 
analysis. The present analysis using running X' is the correct way to reduce renormalisation scale 
dependence in the cross section. We also observe that K factor is much bigger at the Tevatron 
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Figure 2: NLO K-factor and NNLO K-factor are plotted for sneutrino production at the LHC (left 
panel) and the Tevatron Run-ll ( right panel) as a function of its mass. 



compared to that of at the LHC. The reason behind this is attributed to the different behavior of 
parton densities at the Tevatron and the LHC. Note that parton densities rise steeply as x — > and 
fall off very fast as x — > 1, which means the dominant contribution to the production results from 
the phase space region where x ~ t(= m v /S) becomes small. X at Tevatron (0.05 5s T ^> 0.5) is 
larger compared to that at LHC ( 0.007 5s T 5s 0.07) (see also fig.||). Because of this, at Tevatron the 
valence quark initiated processes dominate while gluon and sea quark initiated processes dominate 
at the LHC. As the mass of the sneutrino increases, that is x approaches to unity, the ^-factor at 
Tevatron naturally falls off. At LHC, in the higher mass region (~ 1 TeV), valence quark densities 
start to dominate and hence it stays almost flat compared to Tevatron. We now turn to study the 
impact of the factorisation scale {jip ) and the renormalisation scale (jj, R ) on the production cross 
section. The factorisation scale dependence for both LHC (left panel) and Tevatron (right panel) 
are shown in upper panels of fig. |, for m v = 300 GeV (LHC), m v = 120 GeV (Tevatron). We 
have chosen }JLr = my for both the LHC and the Tevatron. The factorisation scale is varied between 
/if = 0.1 m v and fLp = 10 m v . We find that the factorisation scale dependence decreases in going 
from LO to NLO to NNLO as expected. 

The dependence of the renormalisation scale dependence on the total cross sections for the 
resonant production of sneutrino at the LHC and the Tevatron is shown in the lower panels of fig. 
|3[ Note that the LO is already }Xr dependent due to the coupling X \\Lr)- We have performed this 
analysis for sneutrino mass m v = 300 GeV (LHC), m v = 120 GeV (Tevatron). We have set the 
factorisation scale jip = m v and the renormalisation scale is varied in the range 0.1 < }lR/m v < 
10. We find significant reduction in the /i# scale dependence when higher order QCD corrections 
are included. It is clear from both the panels of fig. || that our present NNLO result makes the 
predictions almost independent of both factorisation and renormalisation scales. 

We could not discuss or show the results of charged slepton due to page limitation. We request 
reader to follow the Ref.[|j"7|]. 

Acknowledgement: Speaker acknowledges the full support of Saha Institute of Nuclear Physics, 
India where this work has been done. Speaker also thanked RADCOR's organiser for inviting him 
to give a seminar in RADCOR 201 1. 
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Figure 3: In the upper panel, sneutrino production cross sections are plotted against the factorisation scale 
jlp with a fixed renormalisation scale pig = m^ for both LHC and Tevatron energies. In the lower panel, they 
are plotted against the renormalisation scale }Xr with a fixed factorisation scale pip = niv for both LHC and 
Tevatron energies. The mass of the sneutrino is taken to be 300 GeV (120 GeV) at LHC (Tevatron). 
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